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Abstract The singularity space-time metric obtained 
by Krori and BaruapQ satisfies the physical require- 
ments of a realistic star. Consequently we explore the 
possibility of applying the Krori and Barua model to de- 
scribe ultra-compact objects like strange stars. For it to 
become a viable model for strange stars, bounds on the 
model parameters have been obtained. Consequences of 
a mathematical description to model strange stars have 
been analyzed. 
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1 Introduction 

In relativistic astrophysics, understanding the nature 
and exact composition of a specific class of compact 
stars which are more compact than ordinary neutron 
stars, has become a field of active research in recent 
years. A neutron star is the final stage of a gravita- 
tionally collapsed star which, after exhausting all its 
thermo-nuclear fuel, gets stabilized by degenerate neu- 
tron pressure. Soon after the discovery of the parti- 
cle 'neutron' by Chadwick, the existence of neutron 
star was predicted. Later on, the concept got observa- 
tional support with the discovery of pulsars [2] . With the 
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progress in our understanding of particle interactions 
at high energy, theoretical modelling of neutron stars 
have improved considerably over the last few decades [3] . 
However, the nature of particle interactions beyond nu- 
clear density is still poorly understood. The conjecture 
that quark matter might be the true ground state of 
hadrons[H[S], has led to the discussions of an entirely 
new class of stellar bodies composed of deconfined it, 
d and s quarks, called 'strange stars'. It is interest- 
ing to note that a strange matter equation of state 
(EOS) seems to explain the observed compactness of 
many astrophysical objects like Her X-l, 4U 1820-30, 
SAX J 1808.4-3658, 4U 1728-34, PSR 0943+10 and 
RX J185635-3754, whose estimated compactness, oth- 
erwise, cannot be explained in terms of a neutron star 
EOS[5||71|5Pl|IuTin]. Though many more exotic phases 
may exist at the interior of such class of stars, in this pa- 
per, we shall restrict our discussions to the strange mat- 
ter EOS only. We shall choose a simple EOS for strange 
quark matter based on the MIT bag model where the 
quark confinement is assumed to be caused by a uni- 
versal pressure B g , called the bag constant 13 . The bag 
model essentially describes the confinement mechanism 
of quarks inside hadrons. By imposing the condition 
that the energy per baryon of strange matter be less 
than that of the nucleon (939 MeV), Farhi and Jaffc[5 
have shown that for a stable strange quark matter the 
bag constant should be B g ~ 60 MeV/fm~ 3 . 

Once the strange matter EOS is known, one can 
employ numerical techniques to get an estimate of the 
gross features of a strange star by integrating the Tolman- 
Oppenheimer-Volkoff(TOV) equations. From the gen- 
eral relativistic view point, question is, what would be 
the relevant back ground space-time to model such class 
of compact stars? In a recent work, Avellar and Horvath[12 
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have considered a wide set of exact and approximate so- 
lutions to model strange stars. The objective of the cur- 
rent investigation is to look for a physically viable ana- 
lytical model which can describe such class of compact 
stars. To this end, we explore the possibility of applying 
the Krori and Barua[T] (henceforth KB) metric to de- 
scribe the interior space-time of a strange star. The KB 
model is singularity free and has earlier been shown to 
be useful to describe realistic starsQi)]. IvanovfTT] has 
shown that the KB model satisfies the necessary en- 
ergy conditions of a realistic star. Varela et aZ[IB], for 
a Einstein-Maxwell system, have used the KB model 
to describe a self-gravitating, charged, anisotropic fluid 
sphere satisfying a linear and/or non-linear EOS. The 
key observations made by them are as follows: (i) spheres 
with vanishing net charge contain fluid elements with 
unbounded proper charge density located at the fluid- 
vacuum interface; (ii) inward-directed fluid forces caused 
by pressure anisotropy may allow equilibrium configu- 
rations with larger net charges and electric field inten- 
sities than those found in studies of charged isotropic 
fluids; (iii) possible applications of the model to de- 
scribe charged strange quark stars, dark matter distri- 
butions and massive charged particles. In a separate pa- 
per, Farook et al[F3\ have used the KB model to analyze 
an anisotropic, charged, static, spherically symmetric 
fluid source. It has been shown that the inclusion of 
a tangential pressure-like variable admits a non-linear, 
Chaplygin-type EOS. Interestingly, the two approaches 
coincide for an EOS of the form p — Hp, where p is 
the energy density, p the pressure and if is a model pa- 
rameter describing the stiffness of the EOS. The results 
obtained by Varela et al may be regained by Farook et 
a/'sQ]5] approach. They too predicted a possible extrap- 
olation of the investigation to the case of astrophysical 
bodies, in particular, for a quark star of radius ~ 8 
km. 

The present investigation is a follow up of the earlier 
works done by Varela et al[TE\ and Rahaman et q/|19]. 
In the present work, we shall incorporate the bag model 
EOS in the KB model and study the subsequent stel- 
lar configurations. It is to be noted here that once we 
specify the EOS, we can integrate the TOV equations 
to derive the gross features of a stellar configuration. 
On the other hand, if one of the metric functions is as- 
sumed a priori, one can determine the subsequent EOS 
of the material composition of the star. However, if both 
the metric functions as well as the EOS are provided, 
it becomes an over determined system. To overcome 
the situation, we include two additional input parame- 
ters into the system. Note that the KB model was orig- 
inally developed for a charged neutral, isotropic and 
spherically symmetric object in static equilibrium. We 



assume here that the composition of the strange star 
is anisotropic in nature coupled with high electric field. 
These assumptions are justified due to the following 
reasons: Strange stars are extremely dense objects and 
at very high density, it is expected that the pressure 
should be anisotropic in nature [Tl]. The electric field 
at the surface of a strange star has also been reported 
to be very high|15l. Implications of these additional pa- 
rameters on the physical behaviour of the strange star 
will be discussed in the following sections. 

The paper has been organized as follows: In Sect. [5J 
we have written the basic field equations. In different 
sub-sections of Sect. [3J we have derived bounds on the 
model parameters based on various physical require- 
ments. In Sect. HI we have discussed implications of 
applying the KB model for the description of strange 
stars. In different sub-sections of Sect. (3J we have ana- 
lyzed various features of the model including stability 
of the resultant configurations. Finally, some conclud- 
ing remarks have been made in Sect. 



2 Basic Equations 

We assume that the interior space-time of a 'strange 
star' is well described by the Krori and BaruajT] metric 
given by 



ds 2 



z v{r Ut 2 + e x{r) dr 2 + r 2 {d6 2 + sin 2 9d(/) 2 ), (I) 



where, A(r) = Ar 2 and v{r) = Br 2 + C. In Eq. (fTJ), A, 
B and C are arbitrary constants which will be fixed on 
the ground of various physical requirements. 

For a static charged fluid source with density p(r), 
radial pressure p r (r), tangential pressure p*(r), proper 
charge density er(r) and electric field E(r), the Einstein- 
Maxwell(EM) equations take the form (we employ the 
geometrized units G = c = 1) 



8irp + E 2 = e 
8irp r - E 2 = e 
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where q(r) is the total charge within a sphere of radius 
r. 

Following the MIT bag model, we take the simple 
form of the strange matter EOS 



1, 

Pr = g(p-4B fl ), 



(6) 
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where, B g is the bag constant. With the choice of the 
above EOS, we have a system of five independent equa- 
tions with five unknown parameters namely, p, p r , pt, 
E(r) and a(r). Substituting the metric potentials given 
by A(r) = Ar 2 , v(r) = Br 2 + C, and their derivatives 
in Eqs. ©-(HO), we obtain 
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The charge density is obtained as 

Ar 2 . 3Ar 2 

cr = — i/j + -^3— [2- (A - W)r 2 ] 
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(11) 
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The charge within a sphere of radius r turns out to be 

(12) 
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The anisotropic stress is obtained as 
1 



A= p t - p r 



2B n 



8nr 2 
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For regularity of the electric field, it must vanish at the 
centre, i.e., 



E 2 (r 
which yields 
B 



0) = - 2 (A- 



B) - 8nB g 
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Substituting the value of B g in Eq. [)14p. we note that 
the parameter A corresponds to the central density 
given by 

87rpo 



.4 



(17) 



Eq. (1171) implies that A is finite and positive. From 
Eq. (Til))) , it then follows that for a positive value of 
the Bag constant, we must have A > B. 

The two pressures and density should be decreas- 
ing functions of r. In our model, radial variation p r is 
obtained as 
dp r 
dr 



B)rAe 



-Ar' 



(18) 



Obviously, at r = 0, ^ = 0. Now, < 0, if the 

condition — A(A+B) < is satisfied. Since A is positive 
and A > B, this implies that B > 0. In Sect. 15.21 
we have shown that if the strong energy condition has 
to be satisfied then B > 0. Therefore, if we assume 
that the matter within the sphere satisfies the strong 
energy condition so that B > 0, then p r will decrease 
radially outward. Similarly, it can be shown that pt also 
decreases radially outward. The radial variation of the 
energy density is obtained as 
dp 3 



■(^4 + B)rAe 



- Ar 



dr 8tt ' 
which also shows that at r 



(19) 



0. 



dp 

dr 



and ^ 



— ^(A + B) < 0. Thus, in this set up, a necessary and 
sufficient condition for regular behaviour of the physical 
parameters will be A > B > 0. 



3 Bounds on the model parameters 



3.2 Regularity at the boundary (r = R): 



One of the advantages of using the KB metric is that 
there there is no singularity in its metric functions. 
Proper bounds, however, should be imposed on the 
conatants appearing in the metric functions so that all 
the physically significant parameters remain well be- 
haved at all interior points of the star. 



3.1 Regularity at the centre (r = 0): 

From Eq. l[7|). we obtain the central density in the form 



Po = p(r = 0) = —(A + B) + B g . 



(14) 



The radius R of the star can be obtained by utilizing 
the condition that the radial pressure should vanish at 
the surface i.e., 



Pr(r = R) = -^{A + B)e 
This yields 



-AR A 



-B g = 0. 



R 




A + B 



(20) 



(21) 



lSirBgi 

Since all the parameters on the right hand side of Eq. (|2"Tj) 
have positive values as discussed in Sect. 13. ll the radius 
of the star is finite and positive. 
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The exterior space-time of the star will be described 
by the Reissner- Nordstrom metric 20 , 21 given by 



ds 2 



2M Q 2 



r r 

2fjn2 i -2 fl jj.2 



+r 2 (d6 z +snT 



dt 2 
), 



2M Q 2 



dr 2 
(22) 



where, Q is the total charge enclosed within the bound- 
ary r — R. Continuity of the metric coefficients gu, g r r 
and across the boundary surface r = R between 
the interior and the exterior regions of the star yields 
the following results: 



1 - 
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Q 2 
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= BRe BR2+c . 


(25) 



Eqs. (|2"3"|) - (f2~5|) determine the values of the constants 
A, B and C in terms of the total mass M, radius R and 
charge Q. By solving the above set of equations, we get 
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(26) 
(27) 
(28) 



Note that the values of the parameters M, R and Q 
should be such that the condition A > B > is satis- 
fied. Moreover, consistency of the above equations puts 
the following constraint on the system: 
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x In 
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= 0. 



(29) 



Eq. (|2"9")l is useful to get an estimate of the charge to 
radius (Q/R) ratio for a given compactness (M/R) of 
the star. Based on logarithmic principle, another con- 
dition that must be fulfilled is that the total charge 
Q 2 < 2RM . Therefore, physical values of the parame- 
ters like the mass, size and charge can not be fixed arbi- 
trarily in this construction. In Sect. 01 we have demon- 
strated that it is possible to choose numerical values of 
masses and radii, consistent with the above constraints. 
These physical values indicate that the corresponding 
configurations are more likely to be strange stars rather 
than neutron stars. 



4 Estimation of physical values 

In this section, we assume the mass and radius of a 
star consistent with the bounds discussed in the previ- 
ous section so that the compactness the star is greater 
than that of a neutron star. This will, in turn, help 
us to get an estimate of the physically relevant pa- 
rameters like the energy density, pressure and the bag 
constant. We have considered compact stars of differ- 
ent compactification factors and calculated the corre- 
sponding constants. The results have been shown in 
Table Q] & [5] For example, for star of mass 1.4 M & 
and radius R = 6.88 km, the values of the constants 
are obtained as A = 0.017977861, B = 0.013506968, 
B g = 0.0002669721 in units of knT 2 and Q 2 / R 2 
[).{)' 11 (Case II). Plugging in G and c in the relevant 
equations, the values turn out to be po — 2.895 x 10 15 
gm cm' 3 , pa — 1.443 x 10 15 gm cm~ 3 , p r (r = 0) = 
p t (r = 0) = 4.361 x 10 35 dyne cm~ 2 and the bag 
constant B g = 202.275 MeV fm~ 3 . Numerical values 
for other cases have been shown in Table [3] Note that 
each case satisfies the condition A > B > 0. Making 
use of the constraint Eq. ([2U)l . we also note that for 
Q 2 /R 2 = 0.004, the minimum value of M/R = 0.25. 
Therefore, for a star of mass 1.4 M , the corresponds 
maximum radius turns out to be R = 8.26 km. To illus- 
trate the behaviour of physical parameters at the inte- 
rior of the star, we have considered case II and plotted 
the variations of the energy density and the two pres- 
sures in Fi. [T| - [31 respectively. 

The bag constant, in this framework, increases with 
the increase of the compactification factor, i.e., the bag 
constant is density dependent. A more dense star re- 
quires a greater bag constant. Similar observations may 
be found in Ref. [35], where a density dependent bag 
constant has been employed to model magnetized strange 
quark stars. The pioneering works of Farhi and Jaffc[5 
showed that for a stable strange matter distribution, 
the bag constant should be around ~ 60 Mev/fm 3 . Our 
results show that with the inclusion of anisotropy and 
electric field, the bag constant turns out to much more 
than its representative value. However, we would like 
to point out here that in Ref. [5], the calculation was 
made for a /3-equilibrium strange quark matter satisfy- 
ing the baryon number conservation principle where the 
charge neutrality condition was employed. The window 
of stability was parametrized by three factors namely, 
the bag constant, the mass of the quark particles and 
the QCD coupling constant. What happens when the 
charge neutrality condition is not imposed is not obvi- 
ous from the analysis. For a relatively massive strange 
quark mass there could be an accumulation of net posi- 
tive charge within the system. Probably, in the presence 



Strange stars in Krori-Barua space-time 



5 



Table 1 Values of Q 2 /R 2 for different choices of the com- 
pactification factor RI/R. 



Case 


M (M Q ) 


R (km) 


M 
R 


Q 1 


I 


1.4 


8.26 


0.25 


0.004 


If 


1.4 


6.88 


0.30 


0.027 


III 


1.4 


5.90 


0.35 


0.061 


IV 


1.4 


5.16 


0.40 


0.105 



Table 2 Values of the model parameters A and B as well as 
the bag constant B g for different choices of compactification 
factor M/R. Data obtained in Case II have been utilized to 
plot the figures. 




Case 


A (km^ 2 ) 


B (km- 2 ) 


B g (km- 2 ) 


I 


0.0102 


0.0073 


0.0001732 


II 


0.01798 


0.01351 


0.000267 


III 


0.0292 


0.0231 


0.0003643 


IV 


0.044 


0.037 


0.000418 



f'r 



Fig. 2 Radial (p r ) and transverse (pt) pressures plotted 
against r. 



Table 3 Energy density, pressure and bag constant for dif- 



ferent 


cases shown in 


Table [T]&[1 








Case 


p(r = 0) 


p(r = R) 


p r (r = 0) 


B a 


0.0020- 




(gm cm" 3 ) 


(gm cm" 3 ) 


(dyne cm -2 ) 


(MeV fm" 3 ) 




I 


1.643 x 10 15 


9.362 x 10 14 


2.123 x 10 35 


131.204 




II 


2.895 x 10 15 


1.443 x 10 15 


4.361 x 10 35 


202.275 


0.0015- 


III 


4.703 x 10 15 


2.015 x 10 15 


8.204 x 10 35 


275.98 


E 1 


IV 


7.087 x 10 15 


2.58 x 10 15 


1.448 x 10 35 


316.699 


0.0010- 



of charge, to counter the repulsive force generated due 
to the electric field, the bag pressure increases. The is- 
sue, however, is a matter of further investigation. What 
we have shown here is that if one wishes to use a mathe- 
matically consistent and physically reasonable analytic 
solution to model strange stars, the bag constant does 
not remain constant. Rather it becomes a free paramter 
which depends on the compactness of the star. 





Fig. 3 The electric field (E 2 ) plotted against r. 

5 Some salient features of the model 

5.1 Generalized TOV equations 

The generalized Tolman-Oppenheimer-Volkoff(TOV) equa- 
tions, in this set up, gets the form 

M G (p + p r ) dp r q X 2 

5 e ^ — - — her — e 2 +-{pt-Pr) = 0,(30) 

r z dr r z r 

where Mq is the effective gravitational mass given by 
M G (r) = \r 2 e^v' = Br z M B ~ A ^- c \ (31) 



Eq. pop describes the equilibrium condition for a charged 
anisotropic fluid subject to gravitational (F g ), hydro- 
static (Fh), electric (F e ) and anisotropic stress (F a ) so 
that 



Fig. 1 The energy density (p) plotted against 



F g + F h + F e + F a = 0, 
where, 

F g = -Br ( p + p r ) = -^-(A + B)e 



(32) 
(33) 
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F h = -^ = -L(A + B)rAe-^ 
ar oir 



(34) 



F P = aEe^r- = 
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27rr 



e- Ar I 2 
2 



A - 3B - — 



(2-{A-3B)i 



xr 2 ( 



-Ar 2 



4irr 2 



(35) 



2 2 

■Fa = - (Pt ~Pr) = ~ 
r r 



1 



8tt 



[(35- 2A + B 2 r 2 - ABr 



2B„ 



(36) 



In Fig. SJ the profiles of these force terms at the in- 
terior of the star for the case II have been shown. The 
plots indicate that an equilibrium stage can be achieved 
under the combined effects of gravitational, electric, hy- 
drostatic and anisotropic stresses. 



Fig. 4 Contributions of different force terms acting on fluid 
elements in static equilibrium. 



5.2 Energy conditions 

The anisotropic charged fluid sphere composed of strange 
matter will satisfy the null energy condition (NEC), 
weak energy condition (WEC) and strong energy con- 
dition (SEC) if the following inequalities hold simulta- 
neously at all points within the star: 

E 2 

" >0, (37) 



8tt 

P + Pr > 0, 

E 2 
E 2 

-p r +2p t + — >0. 

47T 



Pt 



(38) 
(39) 

(40) 



Employing these energy conditions at the centre (r = 
0) , we get the following bounds on the constants A and 
B: 

(i) NEC: p+ f£ >0^> A>0. 

(ii) WEC: p + Pr > A + B > 0, p + p t + g > 
0^A+B>0 

(hi) SEC: p + p r + 2p t + ^ > => B > 0. 
Since the central density is given by po — f^, we must 
have A > 0, i.e., condition (i) is satisfied. The weak and 
strong energy conditions (ii) and (iii) will then be sat- 
isfied if B > 0. With the set of values discussed Sect. |4j 
we have shown in Fig. [5] that the energy conditions are 
satisfied simultaneously within the sphere. 




Fig. 5 Different energy conditions plotted against r for case 
II. 



5.3 Stability 

To examine stability of the resultant configuration, we 
employ the technique based on Herrera's[33] cracking 
(or overturning) concept. Physical acceptability condi- 
tions for a fluid distribution include the condition of 
causality. It suggests that the squares of the radial and 
tangential sound speeds should be within the limit [0, 1]. 
Herrera's [23] cracking (or overturning) concept implies 
that the region for which the radial speed of sound is 
greater than that of transverse speed is a potentially 
stable region. It also suggests that, for 'no cracking' 
to occur, the difference of the two sound speeds, i.e., 
v st ~ v sr should retain the same sign everywhere within 
the matter distribution. In our model, we have 



dp r 

dp 
dp t 
dp 



1 

3' 



-3(A + B)rAe 



-Ar 2 



(41) 
(42) 
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where, 



-At* 



2B 2 r - 2ABr 



(3 = -~2Are 



-Ar J 



B--A + B 2 r 2 - ABr 



1 



2 2 

For causality condition to be satisfied we must have 



< 



-3(j4 + B)rAe~ Ar2 



< 1. 



(43) 



For the assumed set of values, we note that < v 2 sr < 1 
and < v 2 t < 1 as shown in Fig. \6\ Following the 
works of Herrera[53] and Andreasson[26,, we note that 
the configuration will be stable provided | v 2 t — v 2 r \< 
1. Fig. [7]|H] show that the above condition is satisfied 
for the assumed configuration implying stability of the 
configuration. 




Fig. 6 | v^ t — v^ r | plotted against r. 



^0.3,4=0.027| 




Fig. 7 | vf t — fj r | plotted against r. 



Fig. 8 uj t — Vg r plotted against r. 
5.4 Effective mass-radius relation 

For a static spherically symmetric perfect fluid star, 
Buchdahl 24 derived an upper limit for maximum al- 
lowed mass to radius ratio as 2 jf < | which was later 
generalized by Mak et al2h\ for a charged sphere. In this 
model, the effective gravitational mass has the form 

R ( E "\ 2, 1 

p H r ar = - 

o V W 2 

In Fig. [31 we have shown the variation of mass against 
radius. We have also plotted 
which shows that the ratio 
tion of the radial parameter. We note that the con- 
straint on the maximum allowed mass-radius ratio in 
this case is similar to that of an isotropic fluid sphere, 



M, 



eff 



An 



-R{1 



-AR' 



(44) 



M, 



ft 



against R in Fig. [TU1 
M ft f is an increasing func- 



i.e., ^ < |, as obtained by Buchdahl [23]. Defining the 
compactification factor as 

_ M eff (R) 

R 2 

the surface red-shift (Z s ) corresponding to the above 
compactness (u) is obtained as 



AR' 



(45) 



Z s = (1- 2u)-2 - 1 = e~ - 1. (46) 

The maximum surface redshift, in this set up, for a 
strange star of mass 1.4 M Q and radius 6.88 km turns 
out to be Z s = .5303334. 

A lower bound on the mass to radius ratio for a 
charged sphere has been reported by Andreasson[26] 
which has the form 



'M < 



R 



(47) 



IR CP 
3 V 9 3i?' 
This inequality is applicable to stellar objects satisfy- 
ing the inequality p r + 2p t < p. In Fig. 1121 we have 
plotted p r + 2p t — p against r which shows that at 
all interior points, the above condition is satisfied, i.e., 
Andreasson's inequality holds in our model. 



8 



Rahaman et al 




R 

Fig. 9 M e f f plotted against _R. 




Fig. 10 tff plotted against R. 
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Fig. 11 The redshift function Z s plotted against R. 




Fig. 12 2pt + p r — p plotted against r. 



6 Discussions 

We have explored the relevance of KB model pQ in the 
modelling of strange quark stars where the strange mat- 
ter EOS based on the MIT bag model has been as- 
sumed. The inclusion of the EOS does not make the 
system over determined since the matter distribution 
in the set up has been assumed to anisotropic in nature 
together with high electric field. We have shown that a 
self-consistent mathematical model can generate physi- 
cal values which are consistent with a strange star. The 
bag constant, however, in this framework turns out to 
be above the range specified for a stable strange quark 
matter, i.e., (60 — 80 Mev/fm 3 )[SllHl[S]- However, experi- 
mental results from CERN-SPS and RHIC show that a 
wide range of values of the bag constant is possible for a 
density-dependent bag mo del [27]. Perhaps, to compen- 
sate the extra repulsion in the presence of electric field, 
the bag value increases in our set up. The impacts of 
anisotropy and high electric field on the bag model is 
a matter of further investigation and will be taken up 
elsewhere. 
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